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We introduce a time-dependent Eulerian-Lagrangian length-scale and an inverse locality hypoth- 
esis which explain scalings of second order one-particle Lagrangian structure functions observed 
in Kinematic Simulations (KS) of homogeneous isotropic turbulence. Our KS results are consis- 
tent with the physical picture that particle trajectories are more/less autocorrelated if they are 
smoother/rougher as a result of encountering less/more straining stagnation points, thus leading to 
enhanced/reduced turbulent diffusion. 



I. INTRODUCTION 

Turbulent diffusion is a phenomenon of central im- 
portance in oceanic, atmospheric, engineering and as- 
trophysical flows. It is known, since the seminal work 
of Taylor [1], that the turbulent diffusivity (of contam- 
inants, pollutants or other substances advected by the 
turbulent flow) is proportional to the time-integral of 
the Lagrangian velocity autocorrelation function Rl{t) 
of the turbulence. The Lagrangian study of fluid element 
trajectories is therefore central to understanding and cal- 
culating turbulent diffusion. 

A consequence of Taylor's (1921) [1] theory is that 
the turbulent diffusivity is larger/smaller depending 
on whether Lagrangian velocities are correlated over 
longer/shorter times along trajectories. Such long/short 
correlation times reflect large/small Lagrangian veloc- 
ity variations, and Lagrangian velocity variations over 
a time interval r are characterized by the second order 
Lagrangian structure function which is proportional to 
1 — Rl(t). In statistically stationary and isotropic turbu- 
lence it is sufficient to consider only one component of the 
Lagrangian velocity, say v(t); its second order structure 
function is (<5v 2 (t)) = ((v(t) — v(t + t)) 2 ) correspond- 
ing to the Lagrangian velocity power spectrum El(lu) in 
the frequency (w) domain. As noted by Inoue [2] (see 
also Tennekes & Lumley [3]), Kolmogorov similarity ar- 
guments imply that, in the inertial range of times and 
frequencies, (<5v 2 (t)^ ~ er and El(uj) ~ eoj~ 2 where 
e is the kinetic energy dissipation rate per unit mass. 
Support for these scalings has recently been obtained in 
highly turbulent oceanic environments by Lien et al. [4], 
in the laboratory by Mordant et al. [5] and in Direct Nu- 
merical Simulations (DNS) by Yeung [6]. However, Kol- 
mogorov similarity arguments provide very little under- 
standing of the underlying processes responsible for these 
scaling laws. What properties of the spatio-temporal flow 
structure of the turbulence determine that (5v 2 (t)} is 
proportional to r, and what do these properties imply 
for turbulent diffusion? 

In this letter we address this question by use of Kine- 
matic Simulation (KS, see Nicolleau & Vassilicos [7], 



Davila & Vassilicos [8] and references therein). KS is 
a Lagrangian model of turbulent diffusion based on kinc- 
matically simulated turbulent velocity fields which are 
incompressible and consistent with up to second order 
Eulerian statistics of the turbulence such as energy spec- 
tra E(k) in the wavenumber (k) domain. There is no 
assumption of Markovianity or delta correlation in time 
made at any level. Instead, a parameter A controls the 
degree of unsteadiness of the turbulence. It might be 
worth mentioning that, when the energy spectrum input 
has the Kolmogorov fc -5 / 3 form, KS is in good agreement 
with laboratory experiments for one-particle statistics [9] , 
two-particle statistics [7], three-particle statistics [10] and 
concentration variances [11] (particle is used here to mean 
fluid element). KS is also in good agreement with DNS 
for two-particle statistics when the KS energy spectrum 
has the same overall shape as that of the DNS [12]. In 
this letter, we modify the spatio-temporal flow structure 
of the KS turbulence by changing the input energy spec- 
trum E(k) and the input unsteadiness parameter A and 
study how (Sv 2 (t)) changes as a result. 



II. KINEMATIC SIMULATION 

In our KS we use 3D turbulent-like velocity fields of 
the form N k 

u = A„ A k„ cos(k„ • x + uj n t) + 

n=l 

B„ A k„ sin(k„ ■ x + ui n t) (2.1) 

where Nk (typically of order 100) is the number of modes, 
k„ is a random unit vector (k„ = fc n k ra ), and the di- 
rections and orientations of A„ and B„ are chosen ran- 
domly under the constraint that they be normal to k n 
and uncorrelated with the directions and orientations of 
all other wave modes. Note that the velocity field (2.1) 
is incompressible by construction, and also statistically 
stationary, homogeneous and isotropic as shown by Fung 
et al. [13] and Fung & Vassilicos [14]. The amplitudes 
A n and B n of the vectors A„ and B„ are determined 
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by A 2 n = B 2 = |£'(fc„)Afc„ where E(k) is the energy 
spectrum prescribed to be of the form 



E(k) 



v' 2 (p-l) 
2(L/2tt)p~ 1 



(2.2) 



in the range 2tt/L = k\ < k < k^ k = 2ir/r] assumed 
large (i.e. L/rj » 1), and E{k) = otherwise; v' is 
the rms of one velocity component of the KS turbulent- 
like flow; Ak n = (k n+ i — fc„_i)/2 for 2 < n < N). — 1, 
Afci = (fc 2 -fci)/2andAfc A r fc = (k Nk -fcjv fe _i)/2. The dis- 
tribution of wavenumbers is geometric (see [11]), specif- 
ically k n = fcia" -1 with a constant a determined by 
L/i] = a Nk ~ 1 . The frequencies u) n in (2.1) are propor- 
tional to the eddy-turnover frequency of mode n, and 
the dimcnsionless constant of proportionality is A, i.e. 
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FIG. 1. The time evolution of the velocity increment 
(5v 2 (t)^ in a log-log plot. 

Particle trajectories are numerically integrated in (2.1) 
and the velocities of the particles are recorded along their 
flight. By averaging over 5 x 10 3 trajectories in 5 x 10 3 
realizations of the velocity field we obtain (5v 2 (r)) for 
different values of p and A (see Fig. 1 which has been 
obtained for L/rj = 1690, but note that we corroborated 
these results for much higher values of L/rj too). We try 
values of p larger than 1 to ensure that the total kinetic 
energy of the turbulence is finite in the large L/rj limit, 
and values of p smaller than 2 for the power spectrum 
k 2 E(k) of the vorticity and strain rate fields to be an 
increasing function of k (see [8]). 



be t v -C t -C Tl where t v = 2tt /w/v fc - For values of 
A smaller than 1 and close to as well as equal to 0, 
we find (8v 2 (t)) ~ t 2 / 3 over the intermediate range 
t v -C r < T L where T n = 2k / ^ {2it / n)E(2-n: / n) and 
T L = Jq°° Rl(t )dr (t v = t v /X). Fung et al. [13] have 
already observed this 2/3 scaling in frozen (A = 0) KS 
turbulence and added to (2.1) a sweeping of the small- 
scales (2.1) by larger small-scales which restored the Kol- 
mogorov similarity scaling (<5?; 2 (t)) <~ r. Hence, large 
values of the unsteadiness parameter A reproduce this 
sweeping effect on the scaling of (<5v 2 (t)). Kolmogorov 
similarity is conditional on large values of the dimension- 
less parameter A, perhaps a surprising result. 
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FIG. 2. Plot of the change in the power law of the La- 
grangian velocity increment ((5v(t) 2 ) ~ t° as a function 
of p) in KS for different A's. The exponent a and its er- 
ror bars are calculated from a least square fit over the range 
of scales where (5v(t) 2 } ~ t° is observed. This power law 
is harder to detect accurately for values of p closer to 1 be- 
cause a larger proportion of the energy is then distributed to 
higher wavenumbers and both t v and t v increase indefinitely 
as p — > 1 thus requiring an ever larger ratio L/n and therefore 
value of Tl- The fact that more energy is distributed at the 
higher wavenumbers when p decreases means that trajectories 
are more irregular, a fact reflected in the smaller values of a. 
With our extremely good time resolutions and ratios L/rj up 
to 10 4 we still cannot accurately measure values of a below 
about 0.4. 



Running simulations for a variety of values of p be- 
tween 1 and 2 we find 



(Sv 2 (r)) 



~ T 



(3.1) 



III. RESULTS 

A first observation is that when E(k) has the Kol- 
mogorov similarity dependence on wavenumbcr (i.e. p = 
5/3), (<5v 2 (t)) exhibits the Kolmogorov similarity scal- 
ing ^Sv 2 (t)^ ~ r over an intermediate range of time- 
scales only for large enough values of A and not oth- 
erwise. This intermediate range of scales turns out to 



in the intermediate ranges mentioned above and a — 
a(p,X) (see figures 1 and 2). For values of A smaller 
than 1 including A = 0, 

a=p-l, (3.2) 

but for values of A much larger than 1, a is a monotoni- 
cally increasing nonlinear function of p, in fact the same 
function of p for all large enough A, and is such that 
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a > p — 1 and a = 1 at p = 5/3. What is this function 
of pi 

To answer this question we proceed by analogy with 
Richardson's locality hypothesis for two-particle diffusion 
in statistically homogeneous turbulence. According to 
this hypothesis, the rate of change of the mean square 
separation between two particles, i.e. ^ (A 2 ), is a func- 
tion only of (A 2 ) and of the variance of Eulerian velocity 
variations across that length-scale (see [14]). That is 



d_ 

Hi. 



<A 2 ) = / (<A 2 ),(^ 2 (v /<A^)) (3.3) 



where (<5it 2 (r)) is the second order structure function of 
one component of the Eulerian velocity field. It is well 
known that (<5w 2 (r)) <~ r v ~ x for small enough values of r 
is equivalent to E(k) ~ k~ p for large enough values of k. 
The analogy we are trying here, which we call the inverse 
locality hypothesis, is to say that the rate of change of 
((5v 2 (t)) with respect to r is a function only of (<5v 2 (t)) 
and of the length-scale across which the variance of Eu- 
lerian velocity variations is equal to (<5w 2 (r)). That is 



^{5v\r)) = f{{5v\r))Mr)) 



(3.4) 



where £ e (r) is an Eulerian-Lagrangian length-scale de- 
fined as follows: 



(Su 2 (l e (r)))^(Sv 2 (r)). 



(3.5) 



We call l e (j) the equivalent length-scale. Dimensional 
analysis applied to (3.4) and use of (3.5) and (<5w 2 (r)) ~ 
r'P- 1 leads to (3.1) with 



a = 



2(P~1) 
3-p ■ 



(3.6) 



This dependence of a on p fits surprisingly well the results 
obtained by KS for large values of A (see Fig. 2). These 
results reveal, therefore, the existence of an Eulerian- 
Lagrangian length-scale l e (r) which controls the second 
order Lagrangian structure function and Lagrangian fre- 
quency spectrum provided that the sweeping effects are 
taken into account by high values of A. Our results (fig- 
ures 1 and 2 and equations (3.2) and (3.6)) also reveal 
the impact that the two-point spatial flow structure (as 
reflected in the scaling of the energy spectrum) has on 
one-particle Lagrangian statistics whether with or with- 
out sweeping. 

Changes in ((5v 2 (t)) imply changes in Rl{t) and there- 
fore in turbulent diffusion. Taylor's (1921) formula for 
the mean square one-particle displacement (x 2 (t)) is 
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(x 2 (t))=v' 2 [ R L (T)dr (3.7) 
Jo 



which implies a turbulent diffusivity equal to v' 2 T L . Of 
course this relation remains valid in our KS for any value 



of p but the Lagrangian correlation time scale T L is found 
to be an increasing function of p (see Fig. 3 ), thus imply- 
ing a turbulent diffusivity that is also an increasing func- 
tion of p. This conclusion can be explained qualitatively 
by the relation p + 2D s /3 = 3 between p and the frac- 
tal dimension D s of the set of straining stagnation points 
that has recently been obtained by Davila & Vassilicos [8] 
for KS turbulence. As p decreases, D s increases implying 
an increase in the number density of straining stagnation 
points and thereby an increased probability for particle 
trajectories to come in the vicinity of such points and 
experience sudden changes of tack. As a consequence, 
the Lagrangian correlation time scale is smaller/larger 
for smaller/larger values of p. In terms of D s , a = 2 3 ^f s 
when A is large and a = |(3 — D s ) when A is small; a is 
therefore a decreasing function of D s thus reflecting the 
fact that a larger number of changes of tack produces 
more irregular particle trajectories. 
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FIG. 3. Plot of Lagrangian correlation time Tl as a func- 
tion of p for A = 0.5 in KS with L/rj — 100. Similar results 
are found for other values of A and L/rj. 
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